Real-time series formulas for one- and two-compartment pharmacokinetic models

ABSTRACT

The present invention discloses series formulas that can calculate various pharmacokinetic parameters upon administration of multiple maintenance doses and a method to custom-build series equations for dosage regimens that include administration of other doses in addition to the maintenance doses via different routes of administration as a function of dose number and total time in the one- and two-compartment pharmacokinetic models. These series formulas, when compared to the old traditional equations that relate drug concentration to the “dummy” time variable of dosing interval, are ideal for studying time-dependent pharmacokinetic parameters such as the peak time (t max ) after multiple extravascular dosages and they offer an alternative elegant way of writing algorithms for calculating and plotting pharmacokinetic parameters as a function of total time.

BACKGROUND OF THE INVENTION

Mathematical treatments of linear mammillary pharmacokinetic models along with some analytical solutions of drug amount or concentration appeared in the literature in the early 1970 (Wagner, 1981; Benet, 1972, Dubois, 2011). Since then, all mathematical expressions relate to either steady state conditions or time within the first dosing interval (Dipiro et al., 2005). During the last 80+ years there has never been an effort to develop series formulas of drug amount or concentration after administration of multiple doses as a function of real or total time. The present invention relates to a set of mathematical series formulas that can be used to calculate drug amount or concentration upon administration of multiple drug doses of same or different size, via different administration routes as a function of dose number and total time in a one- and two-compartment pharmacokinetic models. The advantage of these real-time series equations as compared to their “dummy” time dosing interval variable counterparts are multiple: First, series equations that allow assessment of pharmacokinetic parameters for dosage regimens that include multiple doses of different size can be custom-produced. Second, time-dependent pharmacokinetic parameters such as peak time after administration of multiple extravascular doses can be simulated as a function of dose number and real time. Third, the area under the curve (AUC) and the average drug concentration can be calculated analytically with definite integration of the series formulas. Fourth, computer programs can be written with nested iterative loops using both dose number and real time (example “program_one_comp_nl_mod”). Fifth, graphs of drug concentration as a function of total time can be plotted without additional lines of code that convert dosing interval into total time.

SUMMARY OF THE INVENTION

The present invention was conceived when the inventor realized that the absorption and elimination rate- and by extension time-dependent variable peak time (t_(max)) after multiple oral doses could not be studied as a function of dose number because all available pharmacokinetic formulas are expressed as a function of dosing interval and not real time. Time-dependent parameters can best be monitored with equations that describe time. In addition, having formulas written in terms of dosing intervals snakes it harder and require additional line of codes to plot drug concentration as a function of time in multiple dosing regimens. The present invention relates to series formulas that can be used to determine drug amount, concentration, area under the curve (AUC) and peak time and peak concentration after administration of multiple drug doses as IV bolus, extravascularly or orally and by intermittent IV infusion (IIV), as a function of real time in a one- and two-compartment pharmacokinetic models. The present invention also describes a method that produces custom-build series equations that can calculate pharmacokinetic parameters as a function of dose number and real time for dosage regimens that include administration of other doses, different in size from the maintenance doses.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 Two-compartment mammillary models

DETAILED DESCRIPTION OF THE INVENTION

The approach for each route of administration and compartment model is as follows:

1. Obtain analytical solutions to differential equations

2. Derive the drug amount or concentration terms of a sequence after administration of multiple doses by applying the principle of superposition

3. Derive the pattern of the sequence

4. Determine the partial sums of the series

5. Derive final series formula of drug amount or concentration as a function of dose number and total time

6. Integrate the series formula to derive AUC and average concentration formulas

7. Differentiate the series formula with respect e to derive peak time (t_(max)) as a function of dose number and use it to obtain series formulas for peak concentration as a function of real time after administration of multiple oral or extravascular drug doses that are characterized by first-order absorption rate kinetics.

8. If multiple doses of different size than the maintenance dose are administered at different times, modify series formulas by reducing the dose number index n by the number of different size doses administered prior to the maintenance doses and by adding terms outside of the series maintenance dose summation formula that account for the drug concentration due to the administration of the other doses.

A. One-Compartment Model I. Intermittent IV Infusion (IIV)

1. Differential equations for inputs and outputs: Input: Drug administered by constant IV infusion (zero-order kinetics Output: First-order kinetics

$\begin{matrix} {{\frac{dx}{dt} = {k_{0} - {k \cdot x}}};{{x(0)} = {{00 \leq t \leq {T\frac{dx}{dt}}} = {{- k} \cdot x}}};{{x(T)} = {{{x\left( {{from}\mspace{14mu}{input}} \right)}T} \leq t \leq \tau}}} & \; \end{matrix}$

Analytical solutions: During Inputs (0≤t≤T): Using the integrating factor method and the initial condition,

${x(t)} = {\frac{k_{0}}{k} \cdot \left( {1 - e^{{- k} \cdot t}} \right)}$ $\begin{matrix} {{C(t)} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot t}} \right)}} & {{{where}\mspace{14mu} C} = {\frac{x}{V}\mspace{14mu}{and}}} & {{CL} = {k \cdot V}} \end{matrix}$

When drug infusion has finished (T≤t≤τ):

x(t)=c·e ^(−k·(t−T)); where c is a constant of integration

${x(T)} = {\frac{k_{0}}{k} \cdot \left( {1 - e^{{- k} \cdot T}} \right)}$ ${k(t)} = {\frac{k_{0}}{k} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({t - T})}}}$ ${C(t)} = {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({t - T})}}} = {{{C(T)} \cdot e^{{- k} \cdot {({t - T})}}} = {C_{\max} \cdot e^{{- k} \cdot {({t - T})}}}}}$

2. Sequence terms for multiple doses administered with a dosing interval τ. Notice that drug concentration during infusion interval is expressed in terms of the minimum concentration from the previous dose whereas drug concentration during the elimination phase is expressed in terms of the maximum drug concentration at the end of the infusion dose. We are also going to call the drug concentration during infusion (0≤t≤T) as C_(a) and during the pure elimination phase (T≤t≤τ) when infusion has stopped as C_(e). The terms of the sequence during the first two infusions are: During the first infusion (n=1):

$\begin{matrix} {C_{a,1} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot t}} \right)}} & {0 \leq t \leq T} \\ {C_{\max,1} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right)}} & {t = T} \\ {C_{e,1} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({t - T})}}}} & {T \leq t \leq \tau} \\ {C_{\min,1} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({\tau - T})}}}} & {t = \tau} \end{matrix}$

During the second infusion (n=2):

$\begin{matrix} {C_{a,2} = {{C_{\min,1} \cdot e^{{- k} \cdot {({t - \tau})}}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot {({t - \tau})}}} \right)}}} & {\tau \leq t \leq {\tau + T}} \\ {C_{\max,2} = {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right)}\left( {1 + e^{{- k} \cdot \tau}} \right)}} & {t = {\tau + T}} \\ {C_{e,2} = {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right)}{\left( {1 + e^{{- k} \cdot \tau}} \right) \cdot e^{{- k} \cdot {({t - \tau - T})}}}}} & {{\tau + T} \leq t \leq {2\tau}} \\ {C_{\min,2} = {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \left( {1 + e^{{- k} \cdot \tau}} \right)}e^{{- k} \cdot {({\tau - T})}}}} & {t = {2\tau}} \end{matrix}$

3. Pattern of Sequence

$\left\{ C_{\max,n} \right\} = \left\{ {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({{({n - 1})} \cdot \tau})}}} \right\}$ $\left\{ C_{\min,n} \right\} = \left\{ {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({{n \cdot \tau} - T})}}} \right\}$ $\left\{ C_{a,n} \right\} = {{\left\{ {{C_{\min,{({n - 1})}} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)}} \right\}\mspace{31mu}{\left( {n - 1} \right) \cdot \tau}} \leq t \leq {{\left( {n - 1} \right) \cdot \tau} + T}}$ $\begin{matrix} {\left\{ C_{e,n} \right\} = \left\{ {C_{\max,n} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau} - T})}}} \right\}} & {{{\left( {n - 1} \right) \cdot \tau} + T} \leq t \leq {n \cdot \tau}} \end{matrix}$

4-5. Partial Sums and final formula of Series

$C_{\max,n} = {{\sum\limits_{n = 1}^{\infty}\;{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({{({n - 1})} \cdot \tau})}}}} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \left( {1 + e^{{- k} \cdot \tau} + \ldots + e^{{- k} \cdot {({{({n - 1})} \cdot \tau})}}} \right)}}$

Multiplying and diving the above by (1−e^(−k·τ)),

$\begin{matrix} {C_{\max,n} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- n} \cdot k \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}} & {{eq}.\mspace{14mu} 1} \end{matrix}$

Using the same approach,

$\begin{matrix} {C_{\min,n} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- n} \cdot k \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({\tau - T})}}}} & {{eq}.\mspace{14mu} 2} \\ {C_{a,n,t} = {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{{- k} \cdot {({n - 1})}}\tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 2})} \cdot \tau} - T})}}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)}}} & {{eq}.\mspace{14mu} 3} \\ {C_{e,n,t} = {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- n} \cdot k \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau} - T})}}}} & {{eq}.\mspace{14mu} 4} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{AUC}_{\max,n} = {{\int_{{({n - 1})} \cdot \tau}^{{{({n - 1})} \cdot \tau} + T}{C_{a}\mspace{11mu}{dt}}} + {\int_{{{({n - 1})} \cdot \tau} + T}^{n \cdot \tau}{C_{e}\mspace{11mu}{dt}}}}} & {{eq}.\mspace{14mu} 5} \\ {{AUC}_{\max,n} = {{T \cdot \frac{k_{0}}{CL}} - {\frac{k_{0}}{k \cdot {CL}} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({{n \cdot \tau} - T})}}}}} & {{eq}.\mspace{14mu} 6} \end{matrix}$

8. Series Formulas for an initial (n=1) Loading IIV dose (D_(L)) followed by multiple IIV maintenance doses (k₀).

For n=1, use the same series formula (eq. 3-4) but replace k₀ with D_(L).

$C_{a,1,t} = {\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot t}} \right)}$ $C_{e,1,t} = {\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot t}} \right) \cdot e^{{- k} \cdot {({t - T})}}}$ $C_{\min,1} = {\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot t}} \right) \cdot e^{{- k} \cdot {({\tau - T})}}}$

For n>1, use the same series formula (eq. 3-4) but add C_(min,1) and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {C_{a,{n > 1},t} = {{{\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({t - T})}}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 2})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 2})} \cdot \tau} - T})}}} + {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)}\mspace{34mu} n}} \in \left\lbrack {2,\infty} \right\rbrack}} & {{eq}.\mspace{14mu} 7} \\ {C_{e,n,t} = {{{\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({t - T})}}} + {{\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau} - T})}}}\mspace{31mu} n}} \in \left\lbrack {1,\infty} \right\rbrack}} & {{eq}.\mspace{14mu} 8} \end{matrix}$

Peak and trough concentrations can be determined from eq. 7 and eq. 8 at t=(n−1)·τ+T and t=n·τ, respectively.

$\begin{matrix} {C_{\max,{n > 1},t} = {{\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({n - 1})} \cdot \tau}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}}} & {{eq}.\mspace{14mu} 9} \\ {C_{\min,n} = {{\frac{D_{L}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot e^{{- k} \cdot {({{n \cdot \cdot \tau} - T})}}} + {\frac{k_{0}}{CL} \cdot \left( {1 - e^{{- k} \cdot T}} \right) \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({\tau - T})}}}}} & {{eq}.\mspace{14mu} 10} \end{matrix}$

An example of a Fortran program is shown below:

Module Shared_data ! Purpose: ! To declare data to share between subroutines implicit none SAVE !Declare parameters integer,parameter:: P=120 ! Infusion time in units of min integer,parameter:: P_h=P/60 ! Infusion time in units of h integer,parameter:: tau_h=12 ! Dosing interval in units of h integer,parameter:: tau=720 ! Dosing interval in units of min real,parameter::DM=596.24 ! Maintenance dose for dosage interval in units of mg real,parameter:: k0= DM/real(P_h) ! Maintenance dose infused in mg/h; DM/P_h real,parameter:: k0_min=k0/60.0 ! Maintenance dose infused in mg/min real,parameter::DL_P=945.00 ! Loading dose for one dosage interval in units of mg real,parameter:: DL=DL_P/real(P_h) ! Loading dose infused in units of mg/h; DM/P_h real,parameter:: DL_min=DL/60.00 ! Loading dose in units of mg/min real,parameter:: CL=2.942 ! Clearance in L/h real,parameter:: CL_min=0.04903 ! Clearance in L/min real,parameter:: k_h=0.07783 ! Elimination rate constant in units h-1 real,parameter:: k=0.001297 ! Elimination rate constant in min-1 integer,parameter::s_t=20 ! loop increment in units of min integer,parameter::istart=1 ! istart of the loop and corresponds to the first infusion integer,parameter::iend=6 ! iend of the loop and corresponds to the last infusion integer,paratneter :incr=1 ! increment of the loop in units of n (infusion #) !Declare variables integer:: n ! Infusion number (times of injection) integer:: t ! Time in minutes real:: Ca ! Drug concentration during infusion (ascending) in mg/L real::Ce ! Drug concentration after infusion is terminated (descending) in mg/L real::Ca_DL ! Drug concentration during infusion (ascending phase) in mg/L with loading dose (DL) real::Ce_DL ! Drug concentration after infusion is terminated with administration of a DL real:: Ca_max=0.0 ! Maximum drug concentration during infusion (ascending) in mg/L real::Ce_min=0.0 ! Minimum drug concentration after infusion is terminated (descending) in mg/L real::Ca_max_DL=0.0 ! Maximum drug concentration during infusion in mg/L with loading dose (DL) real::Ce_min_DL=0.0 ! Minimum drug concentration after infusion is terminated with administration of a DL real::AUC=0.0 ! AUC at the end of dosing interval just prior to the infusion of the next dose real::AUC_DL=0.0 ! AUC at the end of dosing interval just prior to the infusion of the DL in (mg*h)/L real::AUC_ss ! Steady state or max AUC end module Shared_data program one_comp_nl_mod ! Purpose: ! To calculate non-steady state concentrations as a function of time ! in a one-compartment model with and without a loading dose ! The novel long equations were used - manuscript eq. 37-40 except for the first loop ! For n=1 with DL, I have used the equation derived from superposition principle eq. 52 ascending curve only ! Record of revisions: ! Date  Programmer Description of change !

 

! 04/03/20  M. Savva Original Code ! 04/06/20  M. Savva Loops modified to relate directly to infusion time (P), dosing interval (tau) and constant step increment. Nested Loops were adapted to shorten the program and allow increasing n without writing additional code USE Shared_data implicit none  n=1 inf_n_1:  Do t = s_t, P, s:_t  Ca=0.0  Ca=k0/CL*(1-exp(-k*(rea1(t-(n-1)*tau))))+k0/CL*(1-exp(-k*real(P)))&  *exp(-k*(real(t-P-(n-2)*tau)))*(1-exp(-k*real((n-1)*tau)))/(1-exp(-k*real(tau)))  Ca_DL=0.0  Ce_DL=0.0  Ca_DL=Ce_DL*exp(-k*real(t-(n-1)*tau))+DL/CL*(1-exp(-k*real(t-(n-1)*tau)))  write(*,100) ‘n = ’,n,‘t = ’,t,‘C = ’,Ca,‘C_DL = ’,Ca_DL  100 format (A4,I1,5x,A4,I5,2(5x,A4,F8.4))  end Do inf_n_1 elim_n_1:  Do t = P+(n-1)*tau, n*tau, s_t   Ce=0.0  Ce=k0/CL*(1-exp(-k*real(P)))*(1-exp(-k*real(n*tau)))/(1-exp(-k*real(tau)))*exp(-k*real(t-(n-1)*tau-P)))  Ce_DL=0.0  Ce_DL=DL/CL*(1-exp(-k*real(P)))*exp(-k*(real(t-P)))+k0/CL*(1-exp(-k*real(P)))&  /(1-exp(-k*real(tau)))*exp(-k*(real(t-P-(n-1)*tau)))*(1-exp(-k*real((n-1)*tau)))  write(*,101) ‘n = ’,n,‘t = ’,t,‘C = ’,Ce,‘C_DL = ’,Ce_DL  101 format (A4,I1,5x,A4,15,2(5x,A4,F8.4))  end Do elim_n_1 outer: Do n=2,6 inf_n:  Do t = (n-1)*tau, P+(n-1)*tau, s_t  Ca=0.0  Ca=k0/CL*(1-exp(-k*(real(t-(n-1)*tau))))+k0/CL*(1-exp(-k*real(P)))&  *exp(-k*(real(t-P-(n-2)*tau)))*(1-exp(-k*real((n-1)*tau)))/(1-exp(-k*real(tau)))  Ca_DL=0.0  Ca_DL=DL/CL*(1-exp(-k*real(P)))*exp(-k*real(-P+t))+k0/CL*(1-exp(-k*real(P)))/(1-exp(-k*real(tau)))&  *exp(-k*(real(-P+t-(n-2)*tau)))*(1-exp(-k*real((n-2)*tau)))+k0/CL*(1-exp(-k*real(t-(n-1)*tau)))  write(*,102) ‘n = ’,n,‘t = ’,t,‘C = ’,Ca,‘C_DL = ’,Ca_DL  102 format (A4,I1,5x,A4,I5,2(5x,A4,F8.4))  end Do inf_n elim_n:  Do t = +P(n-1)*tau, n*tau, s_t  Ce=0.0  Ce=k0/CL*(1-exp(-k*real(P)))*(1-exp(-k*real(n*tau)))/(1-exp(-k*real(tau)))*exp(-k*(real(t-(n-1)*tau-P)))  Ce_DL=0.0  Ce_DL=DL/CL*(1-exp(-k*real(P)))*exp(-k*(real(t-P)))+k0/CL*(1-exp(-k*real(P)))&  /(1-exp(-k*real(tau)))*exp(-k*(real(t-P-(n-1)*tau)))*(1-exp(-k*real((n-1)*tau)))  write(*,103) ‘n = ’,n,‘t = ’,t,‘C = ’,Ce,‘C_DL = ’,Ce_DL  103 format (A4,I1,5x,A4,I5,2(5x,A4,F8.4))  end Do elim_n end Do outer end program one_comp_nl_d

II. Multiple IV Bolus Doses 1. Differential Equations for Inputs and Outputs:

Input: Drug administered by IV bolus Output: First-order kinetics

$\begin{matrix} {{\frac{dx}{dt} = {{- k} \cdot x}};} & {{x(0)} = D} \end{matrix}$

Analytical solutions:

x(t) = D ⋅ e^(−k ⋅ t) ${C(t)} = {\frac{D}{V} \cdot e^{{- k} \cdot t}}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{n} \right\} = \left\{ {{\frac{D}{V_{d}} \cdot e^{{- k} \cdot t}},{{\frac{D}{V_{d}} \cdot e^{{- k} \cdot t}} + {\frac{D}{V_{d}} \cdot e^{{- k} \cdot {({t - \tau})}}}},\ldots} \right\}$

3. Pattern of Sequence

$\left\{ C_{n} \right\} = \left\{ {\frac{D}{V_{d}} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right\}$

Where, t is the total time.

4-5. Partial Sums and final formula of Series

$C_{\mathfrak{n}} = {{\sum\limits_{n = 1}^{\infty}\;{\frac{D}{V_{d}} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} = {\frac{D}{V_{d}} \cdot e^{{- k} \cdot t} \cdot \left( {1 + e^{k \cdot \tau} + \ldots + e^{k \cdot {({n - 1})} \cdot \tau}} \right)}}$

Multiplying and dividing the above by (1−e^(k·τ)),

$\begin{matrix} \begin{matrix} {C_{n,t} = {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} & {n,{t \in \left\lbrack {0,\infty} \right)}} \end{matrix} & {{eq}.\mspace{14mu} 11} \end{matrix}$

Equations for maximum and minimum drug concentration in the absence of loading doses can be derived from eq. 11 by setting t=(n−1)·τ and n·τ, respectively.

$\begin{matrix} \begin{matrix} {C_{\max,n} = {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}} & {n \in \left\lbrack {0,\infty} \right)} \end{matrix} & {{eq}.\mspace{14mu} 12} \\ \begin{matrix} {C_{\min,n} = {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot \tau}}} & {n \in \left\lbrack {0,\infty} \right)} \end{matrix} & {{eq}.\mspace{14mu} 13} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{AUC}_{n,t} = {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot {\int_{{({n - 1})} \cdot \tau}^{t}{e^{{- k} \cdot {({x - {{({n - 1})} \cdot \tau}})}}{dx}}}}} & \; \\ {{AUC}_{n,t} = {\frac{D}{CL} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot \left( {1 - e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)}} & {{eq}.\mspace{14mu} 14} \\ \begin{matrix} {{{{For}\mspace{14mu} t} = {n \cdot \tau}},{{AUC}_{n,t} = {AUC}_{\max,n}}} & \Rightarrow \end{matrix} & \; \\ {{AUC}_{\max,n} = {\frac{D}{CL} \cdot \left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}} & {{eq}.\mspace{14mu} 15} \\ {C_{{ave},t} = {\frac{1}{\left( {t - {\left( {n - 1} \right) \cdot \tau}} \right)} \cdot \frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot {\int_{{({n - 1})} \cdot \tau}^{t}{e^{{- k} \cdot {({x - {{({n - 1})} \cdot \tau}})}}{dx}}}}} & {{eq}.\mspace{14mu} 16} \\ {C_{{ave},\tau} = \frac{{AUC}_{\max,n}}{\tau}} & {{eq}.\mspace{14mu} 17} \end{matrix}$

8. Series Formulas for an initial (n=1) loading IV bolus (D_(L)) followed by multiple IV maintenance boluses (D). n=1, use the same series formula eq. 11 but replace D with D_(L).

$C_{1} = {{\frac{D_{L}}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} = {\frac{D_{L}}{V_{d}} \cdot e^{{- k} \cdot t}}}$

For n≥2, use the same series formula eq. 11 but add C₁ outside of the summation series formula for the maintenance doses and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {{\begin{matrix} \begin{matrix} {C_{{n > 1},t} = {\frac{D_{L}}{V_{d}} \cdot e^{{- k} \cdot t}}} & {{+ \frac{D}{V_{d}}} \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \end{matrix} & {n \in \left\lbrack {2,\infty} \right)} \end{matrix}t} \in \left\lbrack {0,\infty} \right)} & {{eq}.\mspace{14mu} 18} \end{matrix}$

Equations for maximum and minimum drug concentration in the absence of loading doses can be derived from eq. 18 by setting t=(n−1)·τ and t=n·τ, respectively.

$\begin{matrix} {C_{\max,n} = {{\frac{D_{L}}{V_{d}} \cdot e^{{- k} \cdot {({n - 1})} \cdot \tau}} + {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}}} & {{eq}.\mspace{14mu} 19} \\ {C_{\min,n} = {{\frac{D_{L}}{V_{d}} \cdot e^{{- k} \cdot n \cdot \tau}} + {\frac{D}{V_{d}} \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)} \cdot e^{{- k} \cdot \tau}}}} & {{eq}.\mspace{14mu} 20} \end{matrix}$

III. Multiple Extravascular/Oral Doses

1. Differential equations for inputs and outputs: Input: Drug administered extravascularly with absorption from the site of administration to the central circulation following first-order kinetics Output: First-order elimination kinetics

$\begin{matrix} {{\frac{dx}{dt} = {{k_{a} \cdot x_{v}} - {k \cdot x}}};} & {{x(0)} = 0} \end{matrix}$ $\begin{matrix} {{\frac{{dx}_{v}}{dt} = {{- k_{a}} \cdot x_{v}}};} & {{x_{v}(0)} = {F \cdot D}} \end{matrix}$

Analytical solutions:

$C = {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot t} - e^{{- k_{a}} \cdot t}} \right)}$ $C_{v} = {\frac{F \cdot D}{V_{v}} \cdot e^{{- k_{a}} \cdot t}}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{n} \right\} = \left\{ {{\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot t} - e^{{- k_{a}} \cdot t}} \right)},{{\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot t} - e^{{- k_{a}} \cdot t}} \right)} + {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot {({t - \tau})}} - e^{{- k_{a}} \cdot {({t - \tau})}}} \right)}},\ldots} \right\}$

3. Pattern of Sequence

$\left\{ C_{n} \right\} = \left\{ {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}} - e^{{- k_{a}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)} \right\}$

Where, t is the total time. 4-5. Partial Sums and final formula of Series

$\begin{matrix} {C_{n} = {\sum\limits_{n = 1}^{\infty}\;{\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}} - e^{{- k_{a}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)}}} & \; \\ {{{C_{n,t} = {{\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\{ {{e^{{- k_{a}} \cdot {({t + \tau})}} \cdot \frac{\left( {1 - e^{k_{a} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}} - {e^{{- k} \cdot {({t + \tau})}} \cdot \frac{\left( {1 - e^{k \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}} \right\}}\mspace{31mu} n}},{t \in \left\lbrack {0,\infty} \right)}}{or}} & {{eq}.\mspace{14mu} 21} \\ {C_{n,t} = {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\{ {{e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}} - {e^{{- k_{a}} \cdot {({t - {{({n - 1})} \cdot \tau}})}} \cdot \frac{\left( {1 - e^{{- k_{a}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}}} \right\}}} & {{eq}.\mspace{14mu} 22} \end{matrix}$

The formula for trough drug concentration in the absence of loading doses can be derived from eq. 21 (or eq. 22) at t=n·T.

$\begin{matrix} {C_{\min,n} = {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\{ {{e^{{- k} \cdot \tau} \cdot \frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}} - {e^{{- k_{a}} \cdot \tau} \cdot \frac{\left( {1 - e^{{- k_{a}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}}} \right\}}} & {{eq}.\mspace{14mu} 23} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{AUC}_{\max,n} = {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\lbrack {\frac{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}{k} - \frac{\left( {1 - e^{{- k_{a}} \cdot n \cdot \tau}} \right)}{k_{a}}} \right\rbrack}} & {{eq}.\mspace{14mu} 24} \end{matrix}$

7. Differentiate the series formula with respect to time to derive peak real time (t_(max))and use it to obtain series formulas for peak concentration as a function of real time.

$\begin{matrix} {t = {t_{\max,n} = {{\left( {n - 1} \right) \cdot \tau} + \frac{\ln\left( \frac{k_{a}}{k} \right)}{k_{a} - k} + {\frac{1}{k_{a} - k} \cdot {\ln\left( {\frac{\left( {1 - e^{{- k_{a}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)} \cdot \frac{\left( {1 - e^{{- k} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot n \cdot \tau}} \right)}} \right)}}}}} & {{eq}.\mspace{14mu} 25} \\ {C_{\max,n} = {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\{ {{e^{{- k_{a}} \cdot {({t_{\max,n} + \tau})}} \cdot \frac{\left( {1 - e^{k_{a} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}} - {e^{{- k} \cdot {({t_{\max,n} + \tau})}} \cdot \frac{\left( {1 - e^{k \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}}} \right\}}} & {{eq}.\mspace{14mu} 26} \end{matrix}$

8. Series Formulas for an initial (n=1) loading dose (D_(L)) followed by multiple maintenance doses (D) extravascularly. For n=1, use the same series formula eq. 21 but replace 1) with D_(L).

$C_{1} = {\frac{F \cdot D_{L} \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot t} - e^{{- k_{a}} \cdot t}} \right)}$

For n≥2, use the same series formula eq. 21 but add C₁ outside of the summation series and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {C_{{n > 1},t} = {{{\frac{F \cdot D_{L} \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot t} - e^{{- k_{a}} \cdot t}} \right)} + {{\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left\{ {{e^{{- k} \cdot {({t - {{({n - 1})} \cdot \tau}})}} \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}} - {e^{{- k_{a}} \cdot {({t - {{({n - 1})} \cdot \tau}})}} \cdot \frac{\left( {1 - e^{{- k_{a}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}}} \right\}}\mspace{31mu} n}} \in {\left\lbrack {2,\infty} \right\rbrack\mspace{31mu} t} \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 27} \end{matrix}$

The formulas for peak and trough drug concentration are:

$\begin{matrix} {C_{\min,n} = {{\frac{F \cdot D_{L} \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {e^{{- k} \cdot n \cdot \tau} - e^{{- k_{a}} \cdot n \cdot \tau}} \right)} + {\frac{F \cdot D \cdot k_{a}}{V_{d} \cdot \left( {k_{a} - k} \right)} \cdot \left( {{e^{{- k} \cdot \tau} \cdot \frac{\left( {1 - e^{{- k} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k} \cdot \tau}} \right)}} - {e^{{- k_{a}} \cdot \tau} \cdot \frac{\left( {1 - e^{{- k_{a}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{a}} \cdot \tau}} \right)}}} \right)}}} & {{eq}.\mspace{14mu} 28} \end{matrix}$

B. Two-Compartment Mammillary Pharmacokinetic Model

B.1. Multiple Extravascular/Oral Doses with First-Order Absorption Kinetics and First-Order Elimination from Both Central and Peripheral Compartments (FIG. 1, top left panel) 1. Differential equations for inputs and outputs Input: Drug administered extravascularly/orally with absorption from the site of administration to the central compartment following first-order kinetics Output: First-order elimination kinetics from central and peripheral compartments

$\begin{matrix} {{\frac{dx_{1}}{dt} = {{{- \left( {k_{10} + k_{12}} \right)} \cdot x_{1}} + {k_{21} \cdot x_{2}} + {k_{31} \cdot x_{3}}}};} & {{x_{1}(0)} = 0} \end{matrix}$ $\begin{matrix} {\frac{{dx}_{2}}{dt} = {{k_{12} \cdot x_{1}} - {\left( {k_{20} + k_{21}} \right) \cdot x_{2}}}} & {x_{2}(0)} \end{matrix} = 0$ $\begin{matrix} {{\frac{dx_{3}}{dt} = {{- k_{31}} \cdot x_{3}}};} & {{x_{3}(0)} = {F \cdot D}} \end{matrix}$

Analytical solutions:

${C_{1}(t)} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{1} \cdot e^{{- k_{31}} \cdot t}}} \right)}$ ${A_{1} = \frac{k_{20} + k_{21} - \lambda_{1}}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{2} - \lambda_{1}} \right)}},\begin{matrix} {{B_{1} = \frac{k_{20} + k_{21} - \lambda_{2}}{\left( {k_{31} - \lambda_{2}} \right) \cdot \left( {\lambda_{1} - \lambda_{2}} \right)}},} & {E_{1} = \frac{k_{20} + k_{21} - \lambda_{31}}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{31} - \lambda_{2}} \right)}} \end{matrix}$ λ₁ + λ₂ = k₁₀ + k₁₂ + k₂₀ + k₂₁ λ₁ ⋅ λ₂ = k₁₀ ⋅ k₂₁ + k₁₀ ⋅ k₂₀ + k₁₂ ⋅ k₂₀ ${C_{2}(t)} = {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left( {{A_{2} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{2} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{2} \cdot e^{{- k_{31}} \cdot t}}} \right)}$ ${A_{2} = \frac{1}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{2} - \lambda_{1}} \right)}},\begin{matrix} {{B_{2} = \frac{1}{\left( {k_{31} - \lambda_{2}} \right) \cdot \left( {\lambda_{1} - \lambda_{2}} \right)}},} & {E_{2} = \frac{1}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{31} - \lambda_{2}} \right)}} \end{matrix}$ ${C_{3}(t)} = {\frac{F \cdot D}{V_{3}} \cdot e^{{- k_{31}} \cdot t}}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{1,n} \right\} = \left\{ {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{1} \cdot e^{{- k_{31}} \cdot t}}} \right)},{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{1} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot {({t - \tau})}}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot {({t - \tau})}}} + {E_{1} \cdot e^{{- k_{31}} \cdot {({t - \tau})}}}} \right)}},\ldots} \right\}$

3. Pattern of Sequence

$\left\{ C_{1,n} \right\} = \left\{ {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)} \right\}$

Where, t is the total time. 4-5. Partial Sums and final formula of Series

$\begin{matrix} {C_{1,n} = {\sum\limits_{n = 1}^{\infty}\;{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)}}} & \; \\ {{C_{1,n,t} = {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{31mu} n}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 30} \\ {{C_{2,n,t} = {{\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left\{ {{A_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{2} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{31mu} n}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 31} \end{matrix}$

Trough concentration can be obtained from the general equations above at t=n·τ.

$\begin{matrix} {C_{1,\min,n} = {{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}\mspace{31mu} n} \in \left\lbrack {0,\infty} \right\rbrack}} & {{eq}.\mspace{14mu} 32} \\ {C_{2,\min,n} = {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{1}} \cdot \left\{ {{A_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{2} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}} & {{eq}.\mspace{14mu} 33} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{{AUC}_{\max,n} = {\int_{{({n - 1})} \cdot \tau}^{n \cdot \tau}{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ \ {{A_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}{dt}}}}{{AUC}_{\max,n} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\lbrack {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\lambda_{1}}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\lambda_{2}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{k_{31}}}} \right\rbrack}}} & {{eq}.\mspace{14mu} 34} \end{matrix}$

7. Differentiate the series formula with respect to time to derive peak real time (t_(max)) and use it to obtain series formulas for peak concentration as a function of real time. Due to lack of explicit solutions, the lower eigenvalue exponential term is removed from C_(1,n,t) expression. Given that E₁=−(A₁+B₁) and assuming that λ₁>λ₂ (if λ₁<λ₂ remove λ₁ exponential term from the eq. 30),

$\begin{matrix} {{C_{1,n,t,{apr}} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}}{\frac{{dC}_{1,n,t,{apr}}}{dt} = {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{{- \lambda_{1}} \cdot A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {k_{31} \cdot E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}}} \right\}} = 0}}} & \; \\ {t_{\max,n} = {{\left( {n - 1} \right) \cdot \tau} + \frac{\ln\left( \frac{k_{31} \cdot \left( {A_{1} + B_{1}} \right) \cdot \left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right) \cdot \left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\lambda_{1} \cdot A_{1} \cdot \left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right) \cdot \left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \right.}{k_{31} - \lambda_{1}}}} & {{eq}.\mspace{14mu} 35} \\ {C_{1,\max,n} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}}} \right\}}} & {{eq}.\mspace{14mu} 36} \end{matrix}$

8. Series Formulas for an initial (n=1) loading dose (D_(L)) followed by multiple maintenance doses (D) extravascularly.

For n=1, use the same series formula eq. 30 but replace D with D_(L).

$C_{1} = {\frac{F \cdot D_{L} \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{1} \cdot e^{{- k_{31}} \cdot t}}} \right)}$

For n≥2, use the same series formula eq. 30 (or eq. 31 for the peripheral compartment) but add C₁ outside of the summation formula and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {C_{{n > 1},t} = {{{\frac{F \cdot D_{L} \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{1} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{31mu} n}} \in \left\lbrack {1,\infty} \right\rbrack}} & {{eq}.\mspace{14mu} 37} \\ {C_{2,{n > 2}} = {{\frac{F \cdot D_{L} \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left( {{A_{2} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{2} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{2} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left\{ {{A_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{2} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}}} & {{eq}.\mspace{14mu} 38} \end{matrix}$

Trough concentration can be obtained from the general equations above at t=n·τ.

$\begin{matrix} {C_{1,\min,n} = {{{\frac{F \cdot D_{L} \cdot k_{31}}{V_{1}} \cdot \left( {{A_{1} \cdot e^{{- \lambda_{1}} \cdot n \cdot \tau}} + {B_{1} \cdot e^{{- \lambda_{2}} \cdot n \cdot \tau}} + {E_{1} \cdot e^{{- k_{31}} \cdot n \cdot \tau}}} \right)} + {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{1} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{1} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}\mspace{31mu} n}} \in \left\lbrack {0,\infty} \right\rbrack}} & {{eq}.\mspace{14mu} 39} \\ {C_{2,\min,n} = {{\frac{{F \cdot D \cdot k_{31}}k_{12}}{V_{2}} \cdot \left( {{A_{2} \cdot e^{{- \lambda_{1}} \cdot n \cdot \tau}} + {B_{2} \cdot e^{{- \lambda_{2}} \cdot n \cdot \tau}} + {E_{2} \cdot e^{{- k_{31}} \cdot n \cdot \tau}}} \right)} + {\frac{{F \cdot D \cdot k_{31}}k_{12}}{V_{2}} \cdot \left\{ {{A_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{2} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{2} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}}} & {{eq}.\mspace{14mu} 40} \end{matrix}$

B.2. Multiple IV Bolus Doses with First-Order Elimination from both Central and Peripheral Compartments (FIG. 1, top right panel) 1. Differential equations for inputs and outputs Input: Drug administered by an IV bolus directly into the central compartment Output: First-order elimination kinetics from central and peripheral compartments

$\begin{matrix} {{\frac{dx_{1}}{dt} = {{{- \left( {k_{10} + k_{12}} \right)} \cdot x_{1}} + {k_{21} \cdot x_{2}}}};} & {{x_{1}(0)} = D} \end{matrix}$ $\begin{matrix} {{\frac{dx_{2}}{dt} = {{k_{12} \cdot x_{1}} - {\left( {k_{20} + k_{21}} \right) \cdot x_{2}}}};} & {{x_{2}(0)} = 0} \end{matrix}$

Analytical solutions:

${C_{1}(t)} = {\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)}$ ${A_{3} = \frac{k_{20} + k_{21} - \lambda_{1}}{\left( {\lambda_{2} - \lambda_{1}} \right)}},\mspace{14mu}{B_{3} = \frac{k_{20} + k_{21} - \lambda_{2}}{\left( {\lambda_{1} - \lambda_{2}} \right)}},{{{\lambda_{1} + \lambda_{2}} = {k_{10} + k_{12} + k_{20} + k_{21}}};}$ λ₁ ⋅ λ₂ = k₁₀ ⋅ k₂₁ + k₁₀ ⋅ k₂₀ + k₁₂ ⋅ k₂₀ + k₁₂ ⋅ k₂₁ ${C_{2}(t)} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left( {e^{{- \lambda_{1}} \cdot t} - e^{{- \lambda_{2}} \cdot t}} \right)}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{1,n} \right\} = \left\{ {{\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)},{{\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)} + {\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot {({t - \tau})}}} - {B_{3} \cdot e^{{- \lambda_{2}} \cdot {({t - \tau})}}}} \right)}},\ldots} \right\}$

3. Pattern of Sequence

$\left\{ C_{1,n} \right\} = \left\{ {\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)} \right\}$

Where, t is the total time. 4-5. Partial Sums and final formula of Series

$\begin{matrix} {C_{1,n} = {\sum\limits_{n = 1}^{\infty}\;{\frac{D}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)}}} & \; \\ {C_{1,n,t} = {\frac{D}{V_{1}} \cdot \left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}} & {{eq}.\mspace{14mu} 41} \\ {C_{2,n,t} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left\{ {{\frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {\frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}} & {{eq}.\mspace{14mu} 42} \end{matrix}$

C_(max,n) and C_(min,n) can be obtained from the above formulas at t=(n−1)·τ and t=n·τ, respectively.

$\begin{matrix} {C_{1,\max,n} = {\frac{D}{V_{1}} \cdot \left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}}} \right\}}} & {{eq}.\mspace{14mu} 43} \\ {C_{2,\max,n} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left\{ {\frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} + \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}} \right\}}} & {{eq}.\mspace{14mu} 44} \\ {C_{1,\min,n} = {\frac{D}{V_{1}}\left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}}} \right\}}} & {{eq}.\mspace{14mu} 45} \\ {C_{2,\min,n} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left\{ {{\frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {\frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}}} \right\}}} & {{eq}{.46}} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{{AUC}_{\max,n} = {\int_{{({n - 1})} \cdot \tau}^{n \cdot \tau}{{\frac{D}{V_{1}} \cdot \left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}{dt}}}}{{AUC}_{\max,n} = {\frac{D}{V_{1}} \cdot \left\lbrack {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\lambda_{1}}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\lambda_{2}}}} \right\rbrack}}} & {{eq}.\mspace{14mu} 47} \end{matrix}$

8. Series Formulas for an initial (n=1) loading dose (D_(L)) followed by multiple maintenance doses (D) extravascularly. For n=1, use the same series formula eq. 41 but replace D with D_(L).

$\begin{matrix} {C_{1} = {\frac{D_{L}}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)}} & {{eq}.\mspace{14mu} 48} \end{matrix}$

For n≥2, use the same series formula eq. 41 (or eq. 42 for the peripheral compartment) but add C₁ and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered)

$\begin{matrix} {C_{1,{n > 1},t} = {{{\frac{D_{L}}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{3} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)} + {{\frac{D}{V_{1}} \cdot \left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{14mu} n}} \in \left\lbrack {2,\infty} \right)}} & {{eq}.\mspace{14mu} 49} \\ {C_{2,{n > 1},t} = {{{\frac{D_{L} \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left( {e^{{- \lambda_{1}} \cdot t} + e^{{- \lambda_{2}} \cdot t}} \right)} + {{\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} \cdot \lambda_{1}} \right)} \cdot \left\{ {{\frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {\frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{14mu} n}} \in \left\lbrack {2,\infty} \right)}} & {{eq}.\mspace{14mu} 50} \end{matrix}$

C_(max,n) and C_(min,n) can be obtained from the above formulas at t=(n−1)·τ and t=n·τ, respectively.

$\begin{matrix} {C_{1,\max,n} = {\frac{D_{L}}{V_{1}} \cdot \left( {{A_{3} \cdot e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} + {B_{3} \cdot \left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)} + {\frac{D}{V_{1}} \cdot \left\{ {{A_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)}} + {B_{3} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}}} \right\}}} \right.}} & {{eq}.\mspace{14mu} 51} \\ {C_{2,\max,n} = {{\frac{D_{L} \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left( {e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}} + e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} \right)} + {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left\{ {\frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} + \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}} \right\}}}} & {{eq}.\mspace{14mu} 52} \end{matrix}$

B.3. Multiple Extravascular/Oral Doses with first-Order Absorption Kinetics and First-Order Elimination only from Central Compartment (FIG. 1, bottom left panel)

1. Differential equations for inputs and outputs

Input: Drug administered extravascularly/orally with absorption from the site of administration to the central compartment following first-order kinetics Output: First-order elimination kinetics from central compartment

${\frac{{dx}_{1}}{dt} = {{{- \left( {k_{10} + k_{12}} \right)} \cdot x_{1}} + {k_{21} \cdot x_{2}} + {k_{31} \cdot x_{3}}}};\mspace{14mu}{{x_{1}(0)} = 0}$ ${\frac{{dx}_{2}}{dt} = {{k_{12} \cdot x_{1}} - {k_{21} \cdot x_{2}}}};\mspace{14mu}{{x_{2}(0)} = 0}$ ${\frac{{dx}_{3}}{dt} = {{- k_{31}} \cdot x_{3}}};\mspace{14mu}{{x_{3}(0)} = {F \cdot D}}$

Analytical solutions:

${C_{1}(t)} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{4} \cdot e^{{- k_{31}} \cdot t}}} \right)}$ ${{A_{4} = \frac{k_{21} - \lambda_{1}}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{2} - \lambda_{1}} \right)}};}\mspace{14mu}$ ${{B_{4} = \frac{k_{21} - \lambda_{2}}{\left( {k_{31} - \lambda_{2}} \right) \cdot \left( {\lambda_{1} - \lambda_{2}} \right)}},\;\mspace{11mu}{E_{4} = \frac{k_{21} - k_{31}}{\left( {k_{31} - \lambda_{2}} \right) \cdot \left( {k_{31} - \lambda_{2}} \right)}}}$ λ₁ + λ₂ = k₁₀ + k₁₂ + k₂₁; λ₁ ⋅ λ₂ = k₁₀ ⋅ k₂₁ ${C_{2}(t)} = {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{1}} \cdot \left( {{A_{5} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{5} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{5} \cdot e^{{- k_{31}} \cdot t}}} \right)}$ ${A_{5} = \frac{1}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{2} - \lambda_{1}} \right)}},{B_{5} = \frac{1}{\left( {k_{31} - \lambda_{2}} \right) \cdot \left( {\lambda_{1} - \lambda_{2}} \right)}},{E_{5} = \frac{1}{\left( {k_{31} - \lambda_{1}} \right) \cdot \left( {\lambda_{31} - \lambda_{2}} \right)}},{{C_{3}(t)} = {\frac{F \cdot D}{V_{3}} \cdot e^{{- k_{31}} \cdot t}}}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{1,n} \right\} = \left\{ {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{4} \cdot e^{{- k_{31}} \cdot t}}} \right)},{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{4} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot {({t - \tau})}}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot {({t - \tau})}}} + {E_{4} \cdot e^{{- k_{31}} \cdot {({t - \tau})}}}} \right)}},\ldots}\mspace{11mu} \right\}$

3. Pattern of Sequence

$\left\{ C_{1,n} \right\} = \left\{ {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)} \right\}$

Where, t is the total time. 4-5. Partial Sums and final formula of Series

$\begin{matrix} {{C_{n} = {{\sum\limits_{n = 1}^{\infty}\;{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)}C_{1,n,t}}} = {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{14mu} n}}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 55} \\ {{C_{2,n,t} = {{\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left\{ {{A_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{5} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{14mu} n}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 56} \\ {{C_{3,n,t} = {{\frac{F \cdot D}{V_{3}} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}\mspace{14mu} n}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 57} \end{matrix}$

Trough concentrations can be obtained from the general equations above at t=n·τ

$\begin{matrix} {{C_{1,\min,n} = {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}n}},{t \in \left\lbrack {0,\infty} \right)}} & {{eq}.\mspace{14mu} 58} \\ {C_{2,\min,n} = {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left\{ {{A_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau} \cdot B_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{5} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}} & {{eq}.\mspace{14mu} 59} \end{matrix}$

6. AUC formulas

$\begin{matrix} {{{AUC}_{\max,n} = {\int_{{({n - 1})} \cdot \tau}^{n \cdot \tau}{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}{dt}}}}{{AUC}_{\max,n} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\lbrack {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\lambda_{1}}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\lambda_{2}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{k_{31}}}} \right\rbrack}}} & {{eq}.\mspace{14mu} 60} \end{matrix}$

7. Differentiate the series formula with respect to time to derive peak real time (t_(max)) and use it to obtain series formulas for peak concentration as a function of real time. Due to lack of explicit solutions, the lower eigenvalue exponential term will be removed from expression. Given that E₄=−(A₄+B₄) and assuming that λ₁>λ₂ (if λ₁<λ₂ remove λ₁ exponential term from the equation below),

$\begin{matrix} {{C_{1,n,t,{apr}} = {{{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{31}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\frac{{dC}_{1,n,t,{apr}}}{dt}} = {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{{- \lambda_{1}} \cdot A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {k_{31} \cdot E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}}} \right\}} = 0}}}{t_{\max,n} = {{\left( {n - 1} \right) \cdot \tau} + \frac{\ln\left( \frac{k_{31} \cdot \left( {A_{4} + B_{4}} \right) \cdot \left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right) \cdot \left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\lambda_{1} \cdot A_{4} \cdot \left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right) \cdot \left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \right)}{k_{31} - \lambda_{1}}}}} & {{eq}.\mspace{14mu} 61} \\ {C_{1,\max,n} = {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t_{\max,n} - {{({n - 1})} \cdot \tau}})}}}} \right\}}} & {{eq}.\mspace{14mu} 62} \end{matrix}$

8. Series Formulas for an initial (n=1) loading dose (D_(L)followed by multiple maintenance doses (D) extravascularly. For n=1, use the same series formula eq. 55 but replace D with D_(L).

$C_{1} = {\frac{F \cdot D_{L} \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{4} \cdot e^{{- k_{31}} \cdot t}}} \right)}$

For n≥2, use the same series formula eq. 55 (or eq. 56 for the peripheral compartment) but add C_(L) and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {C_{1,{n > 1},t} = {{{\frac{F \cdot D_{L} \cdot k_{31}}{V_{1}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{4} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {{\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{4} \cdot \frac{\left( {1 - e^{{- k_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{2}} \cdot \tau}} \right)} \cdot e^{{- k_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}n}} \in {\left\lbrack {2,\infty} \right)\quad}}} & {{eq}.\mspace{14mu} 63} \\ {{C_{2,{n > 1},t} = {{\frac{F \cdot D_{L} \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left( {{A_{5} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{5} \cdot e^{{- \lambda_{2}} \cdot t}} + {E_{5} \cdot e^{{- k_{31}} \cdot t}}} \right)} + {{\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left\{ {{A_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{5} \cdot \frac{\left( {1 - e^{{- k_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{2}} \cdot \tau}} \right)} \cdot e^{{- k_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {E_{5} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}n}}},{t \in {\left\lbrack {0,\infty} \right)\quad}}} & {{eq}.\mspace{14mu} 64} \end{matrix}$

Trough concentrations can be obtained from the general equations above at t=n·τ.

$\begin{matrix} {C_{1,\min,n} = {{\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left( {{A_{4} \cdot e^{{- \lambda_{1}} \cdot n \cdot \tau}} + {B_{4} \cdot e^{{- \lambda_{2}} \cdot n \cdot \tau}} + {E_{4} \cdot e^{{- k_{31}} \cdot n \cdot \tau}}} \right)} + {\frac{F \cdot D \cdot k_{31}}{V_{1}} \cdot \left\{ {{A_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{4} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{4} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}}} & {{eq}.\mspace{14mu} 65} \\ {C_{2,\min,n} = {{\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{2}} \cdot \left( {{A_{5} \cdot e^{{- \lambda_{1}} \cdot n \cdot \tau}} + {B_{5} \cdot e^{{- \lambda_{2}} \cdot n \cdot \tau}} + {E_{5} \cdot e^{{- k_{31}} \cdot n \cdot \tau}}} \right)} + {\frac{F \cdot D \cdot k_{31} \cdot k_{12}}{V_{1}} \cdot \left\{ {{A_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{5} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}} + {E_{5} \cdot \frac{\left( {1 - e^{{- k_{31}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- k_{31}} \cdot \tau}} \right)} \cdot e^{{- k_{31}} \cdot \tau}}} \right\}}}} & {{eq}.\mspace{14mu} 66} \end{matrix}$

B.4. Multiple IV Bolus Doses with First-Order Elimination only from Central Compartment (FIG. 1, bottom right panel) 1. Differential equations for inputs and outputs Input: Drug administered by an IV bolus directly into the central compartment Output: First-order elimination kinetics from central compartment

${\frac{{dx}_{1}}{dt} = {{{- \left( {k_{10} + k_{12}} \right)} \cdot x_{1}} + {k_{21} \cdot x_{2}}}};\mspace{14mu}{{x_{1}(0)} = D}$ ${\frac{{dx}_{2}}{dt} = {{k_{12} \cdot x_{1}} - {k_{21} \cdot x_{2}}}};\mspace{14mu}{{x_{2}(0)} = 0}$

Analytical solutions:

${C_{1}(t)} = {\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)}$ ${A_{6} = \frac{k_{21} - \lambda_{1}}{\left( {\lambda_{2} - \lambda_{1}} \right)}},{B_{6} = \frac{\lambda_{2} - k_{21}}{\left( {\lambda_{2} - \lambda_{1}} \right)}},{{\lambda_{1} + \lambda_{2}} = {k_{10} + k_{12} + k_{21}}}$ λ₁ ⋅ λ₂ = k₁₀ ⋅ k₂₁ ${C_{2}(t)} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left( {e^{{- \lambda_{1}} \cdot t} - e^{{- \lambda_{2}} \cdot t}} \right)}$

2. Sequence terms for multiple doses administered with a dosing interval τ.

$\left\{ C_{n} \right\} = \left\{ {{\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)},{{\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)} + {\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot {({t - \tau})}}} - {B_{6} \cdot e^{{- \lambda_{2}} \cdot {({t - \tau})}}}} \right)}},\ldots} \right\}$

3. Pattern of Sequence

$\left\{ C_{1,n} \right\} = \left\{ {\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)} \right\}$

Where, t is the total time. 4-5. Partial Sums and final formula of Series

$\begin{matrix} {\mspace{79mu}{{C_{1,n} = {\sum_{n = 1}^{\infty}{\frac{D}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right)}}}{C_{1,n,t} = {\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot t}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}}}} & {{eq}.\mspace{14mu} 67} \\ {C_{2,n,t} = {\frac{D \cdot k_{12}}{V_{2} \cdot \left( {\lambda_{2} - \lambda_{1}} \right)} \cdot \left\{ {{\frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} - {\frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot t}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}} & {{eq}.\mspace{14mu} 42} \end{matrix}$

C_(max,n) and C_(min,n) can be obtained from the formulas above at t=(n−1)·τ and t=n·τ, respectively.

$\begin{matrix} {\mspace{79mu}{C_{1,\max,n} = {\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}}} \right\}}}} & {{eq}.\mspace{14mu} 68} \\ {C_{1,\min,n} = {\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}}} \right\}}} & {{eq}.\mspace{14mu} 69} \end{matrix}$

The equations for the peripheral compartment C_(2,max,n) and C_(2,min,n) are the same as those derived in section B2 (eq. 44 and eq. 46). 6. AUC formulas

$\begin{matrix} {{{AUC}_{\max,n} = {\int_{{({n - 1})} \cdot \tau}^{n \cdot \tau}{{\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}{dt}}}}\mspace{20mu}{{AUC}_{\max,n} = {\frac{D}{V_{1}} \cdot \left\lbrack {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot n \cdot \tau}} \right)}{\lambda_{1}}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot n \cdot \tau}} \right)}{\lambda_{2}}}} \right\rbrack}}} & {{eq}.\mspace{14mu} 70} \end{matrix}$

8. Series Formulas for an initial (n=1) loading dose (D_(L)) followed by multiple maintenance doses (D) extravascularly. For n=1, use the same series formula eq. 67 (or eq. 42 for the peripheral compartment) but replace D with D_(L).

$C_{1} = {\frac{D_{L}}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)}$

For n≥2, use the same series formula eq. 67 but add C₁ and reduce the index n by the number of loading doses you have administered prior to the maintenance doses (in this case only one D_(L) was administered).

$\begin{matrix} {C_{1,{n > 1},t} = {{{\frac{D_{L}}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot t}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot t}}} \right)} + {{\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot {({t - {({{({n - 1})} \cdot \tau})}}}}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot {({t - {{({n - 1})} \cdot \tau}})}}}} \right\}}\mspace{14mu} n}} \in \left\lbrack {1,\infty} \right)}} & {{eq}.\mspace{14mu} 71} \end{matrix}$

C_(max,n) and C_(min,n) can be obtained from the above formulas at t=(n−1)·τ and t=n·τ, respectively.

$\begin{matrix} {C_{1,\max,n} = {{\frac{D_{L}}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}}} \right)} + {\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)}}} \right\}}}} & {{eq}.\mspace{14mu} 72} \\ {C_{1,\min,n} = {{\frac{D_{L}}{V_{1}} \cdot \left( {{A_{6} \cdot e^{{- \lambda_{1}} \cdot n \cdot \tau}} + {B_{6} \cdot e^{{- \lambda_{2}} \cdot n \cdot \tau}}} \right)} + {\frac{D}{V_{1}} \cdot \left\{ {{A_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{1}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{1}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{1}} \cdot \tau}} + {B_{6} \cdot \frac{\left( {1 - e^{{- \lambda_{2}} \cdot {({n - 1})} \cdot \tau}} \right)}{\left( {1 - e^{{- \lambda_{2}} \cdot \tau}} \right)} \cdot e^{{- \lambda_{2}} \cdot \tau}}} \right\}}}} & {{eq}.\mspace{14mu} 73} \end{matrix}$

The equations for the peripheral compartment are the same as those derived in section B.2 (eq. 50, eq. 52 and eq. 54).

DEFINITION OF TERMS

-   x: drug amount -   τ: dosing interval in units of time -   T: infusion time -   k: first-order elimination rate constant -   k₀: zero-order rate of drug infusion which is also the maintenance     dose (D_(M)) in IIV -   V: volume of distribution -   D: drug dose -   D_(M): drug maintenance dose -   D_(L): drug loading or booster dose -   C_(a,2): Drug concentration and during the second drug infusion -   C_(e,1): Drug concentration when the first drug infusion is finished     during the pure elimination phase -   T≤t≤τ -   x_(v): Drug amount in the extravascular compartment r site of     administration -   k_(a): First-order absorption rate constant -   C_(v): concentration of the extravascular site of administration -   F: bioavailability factor -   k₃₁: First-order absorption rate constant in the two-compartment     models

REFERENCES

-   Benet, L. Z., General treatment of linear mammillary models with     elimination from any compartment as used in pharmacokinetics. Pharm.     Sci., 61 (4), 536-541,1981. -   J. T., Spruill, W. J, Wade, W. E., Blouin, R. A. and Pruemer, J. M.,     Concepts in Clinical Pharmacokinetics. 2005, American Society of     Health-System Pharmacist, MD, USA. -   Wagner G J, History of Pharmacokinetics, Pharmac. Ther. 12,     537/-562,1981 -   Dubois, A., Mathematical expressions of the pharmacokinetic and     pharmacodynamic models implemented in the PFIM software. 2011.     http:/www.pfim.biostat.fr/PFIM PKPD library.pdf 

I claim:
 1. The real time series formulas in a one-compartment and two-compartment pharmacokinetic models developed and numbered in this document that can be used to: execute an algorithm programmed in high- or low-level language to calculate drug amount, concentration, minimum and maximum drug amounts and concentrations, AUC, average drug concentration as the AUC divided by the dosing interval, peak time and peak concentration; execute an algorithm programmed in a computer language and plot the relationship between the aforementioned calculated pharmacokinetic parameters and as a function of time and dose number. produce series formulas with rearranged numerators and denominators but are otherwise the same formulas. For example, eq. 21 is the same as eq.
 22. 2. A method to custom-build series formulas for administration of other doses that differ in size from the maintenance doses, to calculate all pharmacokinetic parameters described in claim 1 by: adding the concentration(s) or amounts due to administration of the other doses using the principle of superposition outside of the summation series formulas and subtracting the number of other doses from the index n of the original maintenance dose series formulas. 